Ideal cylindrical cloak: Perfect but sensitive to tiny perturbations 
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Abstract 

A cylindrical wave expansion method is developed to obtain the scattering field for an ideal 
two-dimensional cylindrical invisibility cloak. A near-ideal model of the invisibility cloak is set up 
to solve the boundary problem at the inner boundary of the cloak shell. We confirm that a cloak 
with the ideal material parameters is a perfect invisibility cloak by systematically studying the 
change of the scattering coefficients from the near-ideal case to the ideal one. However, due to the 
slow convergence of the zero*' 1 order scattering coefficients, a tiny perturbation on the cloak would 
induce a noticeable field scattering and penetration. 
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The exciting issue of exotic materials invisible to electromagnetic (EM) waves was dis 



cussed in recent works 
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111 ]. Based on a coordinate transformation 



of Maxwell's equations, Pendry et al. first proposed an invisibility cloak, which can protect 
objects inside the cloak from detection When EM waves pass through the invisibility 
cloak, the cloak will deflect the waves, guide them around the object, and return them to 
the original propagation direction without perturbing the exterior field. Numerical methods 
have been applied to solve the EM problem involving invisibility cloaks 6|, |9j, and an experi- 
mental result of the invisibility cloak using metamaterial with simplified material parameters 
has also recently been reported Q. Yet, the ideal invisibility cloak has not been confirmed 
as a perfect cloak, due to the extreme material parameters required (zero or infinity) in the 
ideal cloak when approaching the inner boundary. Also, numerical methods usually describe 
the material parameters discretely, which can be computationally intensive in extreme cases. 
Thus it is preferable to use an analytical or semi-analytical method whenever possible. 

In this paper, we will study the scattering for an ideal invisibility cloak. We focus 
our analysis on the 2D cylindrical cloak, because the wave equation can be simplified in 
comparison with the 3D case, and a 2D invisibility cloak is more feasible to fabricate Q. 
Here we take advantage of the cylindrical geometry of the structure and use the cylindrical 
wave expansion method to study the device semi-analytically. To avoid extreme values 
(zeros or infinity) of material parameters at the cloak's inner surface, we introduce a small 
perturbation into the ideal cloak, and allow the perturbation to approach zero to study the 
scattering problem for the ideal cloak. Such an asymptotic analysis not only can confirm 
whether the ideal cloak would be perfectly invisible or not, it also provides hints on how 
sensitive such a device is to finite perturbations. A sensitivity analysis of the invisibility cloak 
directly determines the possibility of its application. Our studies show that the cylindrical 
invisibility cloak is very sensitive to tiny perturbations of the material parameters. 

First, let's look at the wave equation inside a cylindrical cloak. According to Ref. [lj], a 
simple transformation 



r 



I b— a 



+ a, 9' = 6, z' = z (1) 



can compress space from the cylindrical region < r < b into the annular region a < r' < b, 
where a is the inner radius of the cloak, b is the outer radius of the cloak, and r, 9 and z 
(r', 9' and z') are the radial, angular and vertical coordinates in the original (transformed) 
system, respectively. Following the approach in Ref. jl|, the permittivity and permeability 



tensor components for the cloak shell can be given as 

e r = fi r = J^a, e e = fi e = 

t z f^z \b-a) r ' 

and air is assumed for the ambient environment and the interior regions. In the following, 
the transverse-electric (TE) polarized electromagnetic field is considered (i.e. the electrical 
field only exists in the ^-direction) , however the transverse-magnetic derivation follows in 
similar manner. Throughout the paper, a exp(— iut) time dependence is assumed. For the 
TE-polarized wave, only e z , fi r , and fig are relevant to the following general wave equation 
governing the E z field in the cloak's cylindrical coordinate 

1 8 r 8E 18 1 9E 

e z rdr [ fi e dr } + e z r 2 86 [ fi r 89 } + ^ ~ U ' {6) 

where ko is the wave vector of light in vacuum. If we substitute Eq. [2] for e z , fi r , and fig, we 
find 

Equation |4ban be solved by a separation of variables E z = ty(r)Q(6) and the introduction 
of a constant /: 

(r - af^ + (r - a)^- + [(^_) 2 (r - a) 2 A; 2 - = (5) 

|^ + ' 2 e = o. (6) 

Equation [5] is the Z fh -order Bessel differential equation, and the general solution of Eq. [6] is 
exp(il9). Therefore, there exists a simple set of solutions to E z in the cloak shell of the form 

F,(fci(r - a)) exp(iZ0) (7) 

where fci = kob/(b — a), is the /-order Bessel function, and / is an integer number as 
required by the rotational boundary condition. 

Let us consider the scattering problem in which an arbitrary wave is incident on the 
cloak. According to the rigorous scattering theory [3j , the incident field in the 2D case can 
be expanded in the cloak's coordinates with the following expression 

^ n = X>|V,(A:or)exp(iZ0), (8) 
i 



where Ji is the / -order Bessel function of the first kind. The scattering field can also be 
expanded as 

^r = ^«r^(v)ex P (^), (9) 
i 

where Hi is the Z ift -order Hankel function of the first kind. 




FIG. 1: The schematic of a near-ideal invisibility cloak: The distribution of the material parameter 
is the same as the ideal one shown in Eq. [21 and the outer boundary is still fixed at r = b. However, 
the actual inner boundary is at r = a + 5, where 5 is a very small positive number. 

We note that the scattering coefficients cannot be directly obtained for the ideal cloak 
since e z — > 0, fi r — > 0, and fig — > oo when r — > a, and the Bessel function of the second 
kind in Eq. has a singularity at r = a. In order to circumvent this, we introduce a small 
perturbation to the ideal cloak which we refer to as the near-ideal cloak, see Fig. [TJ We 
expand the inner boundary of the cloak shell slightly, so that it is located at r = a + 5, where 
5 is a very small positive number. However, the material parameters are still calculated 
according to Eq. [2] as if the inner boundary is unchanged. The outer boundary remains 
fixed at r = b. When 5 — > 0, our model will be equivalent to the ideal cloak. Now the 
electric-field in each region can be given by 

(6 < r) E z = Yl a\ n Ji{k a r) exp(i70) + a^ c Hi(k r) exp(iW) 

i 

(a + 5 <r <b) E z = ]T a]Ji(ki(r - a)) exp(z/0) + afH^k^r - a)) exp(z/0) (10) 

i 

(r < a + 5) E z = £ afJi(k r) exp(U9) 

i 

where a\(i = 1,2, 3) are the expansion coefficients for the resulting field inside the cloak. 
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The tangential fields E z and Hg (which can be obtained from E z ), should be continuous 
across the interfaces at r = a + 5 and r — b; and the orthogonality of exp(U8) allows waves 
in each Bessel order to decouple. Thus, we can have the following four equations: 



a^Mkob) + of #,(fc &) = ajj^hib - a)) + afH^b - a)) 
ajJiihd) + afHiihS) = afJi(k (a + 5)) 



(11a) 



fcoof J[{k Q b) + koa^H'^kob) 
ki 



ki 



(lib) 

ajJ'^hib - a)) + -^-afHKhib - a)) (11c) 



ki 



-afH'^k^) = koafJ^koia + 5)) 



;iid) 



He{a + 5) ne{a + 5) 

which is a set of linear equations. Thus each order expansion coefficient in each material 
region can be exactly solved. In turn we can obtain the fields in each region. 



As a direct result of this set linear equations, we can prove that when 5 — > 0, af 



O; 







a; 



a] n , and af — > for any af 1 , i.e., the ideal cloak is a perfect invisibility cloak. Firstly, it 
can be assumed that \a\(i = sc, 1, 2, 3)| must be finite. Otherwise, the scattering field would 
be infinite if the incident field has the I th order component. Secondly, due to ki(b — a) = kob 
and kl = kofig(b), when 5 — > 0, Eq. lllal and lllcl become (a™— a}) Ji(kob) + (ai c —af)Hi(kob) = 
and (af 1 — a})J[{kob) + (af c — a^H'^kob) = 0, respectively. Since b can be arbitrary and 
the Bessel functions are not always zeros, a\ n = a} and af c = af must be satisfied. Thirdly, 
from Eq. Illbl we can obtain the following inequality 



U 2 #,(M)| < \afJi(k (a + 5))\ + la, 1 Ji(h8) I . 



(12) 



When 5^0, the right side of the above inequality approaches a finite value but \Hi(ki5)\ 
approaches infinity on the left side. Thus, \af\ must approach zero. Finally, from Eq. Illbl 
we can also obtain that lafH'^hS)] < af^Ji(k (a + 5)) +\ajJ' l (k 1 S)\. While from Eq. Hid] 
we have 



\k af Ji(k (a + 8))\ < 



of JIM 



Hg(a + 5) 



afHKhS) 



(13) 



Since fig (a + 5) — » oo and the right side of the above inequality approaches zero when 5^0, 
we obtain that \af\ — > 0. Consequently, this argument proves that the scattering field and 
the field in the interior region of the cloak are zero when 5 = 0, i.e., the ideal cloak is a 
perfect invisibility cloak. 
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FIG. 2: (Color online) (a) Snapshot of the resulting electric-field distribution, (b) the corresponding 
norm in the vicinity of the cloaked object, and (c) the snapshot of the corresponding scattering field 
outside the cloak for the near-ideal cloak with 5 = 10~ 5 a and when a plane wave is perpendicularly 
incident on the cloak. The black lines outline the cloak shell. Axis unit: meter. 

Although we have just confirmed that the ideal cloak can provide perfect invisibility, 
further study the near-ideal cloak by the above analytical method illuminates how sensitive 
the parameter 5 is to the performance of the cloak. As an example, we use the same material 
parameters in Ref. jg] where the inner radius of the cloak is a = 0.1m, the outer radius of 
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the cloak is b = 0.2m, and the frequency of the incident plane wave is 2GHz. Similarly, we 
also consider a plane wave incident on the cloak, where the expansion coefficients in Eq. [8] 
are 

a z m = i l Aexp(-ik r l cos(tp + 9i) - Hip), (14) 

where (ri,0i,O) is the coordinate of the phase reference point, A is the amplitude of the 
plane wave, and tp is the incident angle 13J. Here the phase reference point is set at T\ = 4a 
and 6\ = it, the amplitude is A = 1, and the incident angle is tp = (i.e. the plane wave 
propagates from left to right). We use 31 Bessel terms ( — 15 < I < 15) to calculate the 
scattering field for the near-ideal cloak with 5 = 10~ 5 a. The number of expansion terms 
is sufficient for convergence of the calculated fields. Figure [2] shows the snapshot of the 
resulting electric-field distribution (i.e. the real part of the electric-field phasor), and the 
corresponding norm in the vicinity of the cloaked object. The electric-field distribution 
clearly demonstrates the cloaking effect of the near-ideal cloak to the incident plane wave. 
However, the norm of the electric- field (Fig. [2] (b)) reveals that there is still a little bit of the 
field in the cloak interior and an obvious scattering ripple around the cloak. The amplitude 
of the resulting electric-field at the center is 0.197. The snapshot of the scattering field 
(Fig. [2] (c)) shows that it propagates almost isotropically in all angles. Even though the 
amplitude of the scattering field is much smaller than that of the incident plane wave, the 
interference of the incident plane wave and the scattering field creates the ripples in the 
norm (Fig. 12(b)). 

Since each order expansion coefficient of the scattering field is only relevant to each order 
expansion coefficient of the incident field (cf. Eq. [TTT) . we can define the scattering coefficient 
for each order as 

c ? = (15) 

These coefficients for the field inside the cloak c z = a|/a| n , i = 1,2,3 can also be defined 
in the same way. To study the ideal cloak, we more 5 closer to 0. The amplitude and the 
phase of these coefficients for 10~ 8 a < 5 < 10~ 2 a are shown in Fig. [3], where (a)-(b) and 
(c)-(d) correspond to the cases of I = and 1 = 1, respectively. 

From Fig. [31 it is clear that c| is always equal to 1 for both cases. That is, the incident 
field propagates into the cloak without any reflection at the outer boundary, which coincides 
with the explanation of the cloaking effect from the coordination transformation approach 
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FIG. 3: (Color online) The amplitude and phase of the scattering coefficients for the different 5, 
where (a)-(b) and (c)-(d) correspond to the cases of I = and / = 1 respectively. cf c is denoted 
by the blue point-dashed line, Cj (the solid black line), cf^ (the red circle-marked) , and c| 3 ^ (the 
green star-marked). 



in 



1|. The same behavior for the scattering fields occurs at the outer boundary, where they 
propagate from inside to outside without any reflection, thus c\ c is always equal to c\ . 



Our computational results also confirm that both cf c and c\ approach zero when 5 — ► 0. 



(3) 

I 

(3) 



In particular, compared with the case of / = 0, |cf c | and cf' for I = 1 are much smaller, 
and approach zero more rapidly. This is also observed for the other higher order cases. 
Thus, in the case of the plane wave incident, where \a\ n \ is the same for each order, the 
dominating term of the scattering field outside of the cloak is of the form of the zero^-order 
Hankel function of the first kind. Meanwhile, the resulting field in the interior region has a 
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dominating term of Jo(fco^) • This explains the near azimuthally invariable distribution of 
the field in the interior region (see Fig. [21(a)) and the scattering field outside the cloak (see 
Fig. [21(c)), which has also been mentioned in Ref. 6j. 

It is worth noting that the zero*' 1 order scattering coefficients Cq c and decrease ex- 
tremely slowly with reduced 5, e.g. when 5 is decreased from 10~ 5 a to 10 _8 a, |cq c | decreased 
only from 0.175 to 0.099. By utilizing the arbitrary calculation precision of the software 
MATHEMATICA, we found that the convergence of the limit is so slow that even for 
5 = 10~"a (i.e. e z ~ 4 x 10~", /i r ~ 10~", and fig ~ 10" at the inner boundary in this 
case), |cq c | = 6.973 x 10~ 3 . Therefore, we conclude that even though an cloak with the 
ideal material parameters in Ref. jl| is a perfect cloak, a non-ideal invisibility cloak does not 
provide a good enough cloaking effect due to the slow convergence of |cq c | and . 

In conclusion, we have used the cylindrical wave expansion method to study the electro- 
magnetic scattering properties of a 2D invisibility cloak. A near-ideal model of the invisibility 
cloak is set up to solve the boundary problem at the inner boundary of the cloak shell. By 
systematically studying the change of the scattering coefficients from the near-ideal case 
to the ideal one, we confirm that the cloak with the ideal material parameter is a perfect 
invisibility cloak. But due to the slow convergence of the scattering coefficients, a tiny per- 
turbation on the cloak would induce a noticeable field scattering and penetration. We also 
proved that the scattered and penetrated fields are dominated by zero^-order cylindrical 
waves. Though our work has focused on the 2D cylindrical cloak, it can be reliably ex- 
tended to the 3D spherical case. Our method and results are also useful for either designing 
or detecting this type of the invisibility cloak. 
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